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Abstract 

We introduce a new set of flow parameters to describe the time dependence of the 
equation of state and the speed of sound in single field cosmological models. A scale 
invariant power spectrum is produced if these flow parameters satisfy specific dynamical 
equations. We analyze the flow of these parameters and find four types of fixed points that 
encompass all known single field models. Moreover, near each fixed point we uncover new 
models where the scale invariance of the power spectrum relies on having simultaneously 
time varying speed of sound and equation of state. We describe several distinctive new 
models and discuss constraints from strong coupling and superluminality. 
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1 Introduction 



Our Universe is observed to be homogeneous and isotropic at large distances from super- 
galactic scales (~ 1 Mpc) to the whole visible Universe (~ 10 4 Mpc). Moreover, the power 
spectrum of the primordial density fluctuations on such large scales is measured to be nearly 
scale invariant. According to the current view, these perturbations arose as quantum fluc- 
tuations in the early Universe, then left the horizon before the standard expansion phase, 
only to recently reenter. The evolution of the adiabatic fluctuations in the early Universe are 
determined by the equation of state w and the speed of sound c s . The question is, which 
types of early cosmic evolution w(t), c s (t) would give rise to the observed 10 e- folds of nearly 
scale invariant perturbations? 

In a common class of models for generating scale invariant perturbations, the Universe is 
dominated by a scalar field with the canonical kinetic term and certain potential term. In this 
case c s equals 1, and the evolution of the scale factor a, hence H = a/a and e = —H/H 2 = 
|(1 + w), is determined by the form of the potential. Two well known such models are 
the inflation and the adiabatic ekpyrosis, which can both produce a scale invariant power 
spectrum. In the inflationary scenario [HO [3], e stays nearly constant and close to 0, while 
the scale factor a grows exponentially, pushing a wide range of length scales outside the 
horizon. In the adiabatic ekpyrotic scenario [H [5], however, the scale factor a contracts 
slowly while e increases rapidly, so that the horizon scale shrinks and leaves a wide range of 
perturbation modes outside. The scale invariant modes created in a contraction phase can 
be carried on to the standard expansion phase through a cosmic bounce 0, El [H] . 

For a constant c s , the scale invariance of the power spectrum fully constrains the time 
dependence of the scale factor a(t) [9], [10] . In addition to inflation and adiabatic ekpyrosis, 
there is a third type of solution in which the universe goes through an "apex" [9], a slow 
transition from expansion to contraction. Similar to adiabatic ekpyrosis, this scenario also 
depends on a rapidly changing e while the scale factor a is nearly constant. One difference 
is that here the scale invariant modes are generated during the slow expansion phase before 
the apex, instead of during a slow contraction phase. 

More generally, the speed of sound c s may vary with time, as in single field models with a 
noncanonical kinetic term. Several models involving a time varying c s have been considered. 
In the tachyacoustic expansion scenario (111 \12\ [T3l I14j . e is assumed to be constant, and 
the scale invariant power spectrum is produced by having a suitable time dependence of the 
speed of sound c s . Another example is considered in |15| where c s is taken to be proportional 
to e, so that the three-point function of the curvature perturbation is also scale invariant. 

With a time dependent c s , the requirement of a scale invariant power spectrum alone does 
not fully determine the cosmic evolution. There exist infinitely many combinations of a(t) 
and c s (t) that can all produce scale invariant perturbations. We would like to know if these 
different scenarios can be described and categorized into distinctive classes of cosmological 
models. 

In this paper we present a general scheme to describe single field cosmological models 
by a new set of flow parameters. Any single field model can be represented by a point or 
a trajectory in this parameter space. The requirement of a scale invariant power spectrum 
constrains the dynamics of the flow parameters and determines particular classes of trajec- 
tories. We show that all existing models correspond to nearly constant values of the flow 
parameters. By looking for fixed points in the parameter space, we find four general types 
of models that encompass inflation, adiabatic ekpyrosis, apex, and tachyacoustic expansion 
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respectively. Furthermore, we analyze the flow lines near the fixed points in different cross 
sections of the parameter space. By looking for various behavior of each fixed point, we 
obtain new cosmological models that can produce scale invariant perturbations equally well. 

Our approach is very different from the conventional way of starting from a particular 
scalar field Lagrangian and studing its cosmological solutions. In that approach, generally 
speaking, only solutions in a limited region of field phase space would yield scale invariant 
perturbations, corresponding to a particular type of mechanism represented by some part of 
the flow parameter space. In our approach, however, by tracing the flow lines in the whole 
parameter space, we obtain all possible cosmic evolutions that can create a scale invariant 
power spectrum. Each and every flow line leads to scale invariant perturbations, but different 
trajectories may correspond to very different scalar field models. This approach effectively 
avoids the limitations of particular scalar field models. 

In section [2] we define the flow parameters and specify their dynamics in order to produce 
scale invariant perturbations. We use these parameters to study the c s = 1 case and compare 
our results with known models. In section [3] we generalize to a time dependent speed of 
sound and analyze all four fixed points of the flow parameters. In particular, we uncover new 
scenarios and describe the corresponding new types of cosmological models. A few examples 
are shown in section |4j Section [5] provides further discussions on physical constraints from 
strong coupling and superluminality etc. We conclude with a summary in section [6} 

2 Flow parameters 

Consider a single scalar degree of freedom evolving on a Friedmann- Robertson- Walker back- 
ground. The adiabatic fluctuations can be studied by using the effective action for the 
curvature perturbation ( [16\ HOj. At quadratic order, the action is given by 



where q 2 = a 2 e/c s describes an effective background geometry, and ' denotes the derivative 
with respect to the "sound horizon time" y, related to the physical time t by dy = (c s /d)dt. 
Normalizing the field £ by v = V2qC, the quadratic action becomes 



The variable v can be canonically quantized in the standard way [TTj. The Fourier modes 
vt should satisfy the equation of motion 



In comparison to the comoving scale 1/k, the inverse square root of \q"/q\ represents the 
freeze-out horizon size. In general, if q oc y n , then q"/q = n(n — l)/y 2 , hence the horizon size 
is roughly ~ \y\. Modes deep inside the horizon, with k\y\ ^> 1, are effectively in a Minkowski 
background. Specifically, eq. (pi) has an exact solution in terms of the Hankel function, 




(1) 






(3) 




(4) 
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where v = | n — ^ | . This solution has been chosen to match the Minkowski vacuum state in 



the asymptotic past when the mode is deep inside the horizon, Vk ~ e~ lky /\/2k as ky — > — oo. 
The horizon crossing happens at a later time when k\y\ ~ 1. Accordingly, the time variable 
y runs from — oo to 0. 

At late times when k\y\ <C 1, the curvature perturbation = Vk/V%q becomes 



C fc = CiAT^-jOs-"-" + C 2 kr(-y)*- n +" + 0{k 2+ »), (5) 

where C\,C2 are dimensionless Taylor coefficients. Note that for n < ^, the leading term in 
eq. ([5]) is time independent while the second and higher terms are decaying at late times, hence 
Ck approaches constant values on large scales outside the horizon. But for n > ^, however, 
the leading term keeps growing even after crossing the horizon, indicating an instability in 
the curvature perturbation against the homogeneous background. 
The power spectrum of the curvature perturbation £ is given by 

^c = ^iai 2 , (6) 

where Cfc should be evaluated at horizon crossing or when the super-horizon modes stop 
growing [16]. In either cases, the spectral tilt is given by the leading term in ([5]), 

n s - 1 = 3- 2u = 3- |2n- 1|. (7) 

For the power spectrum to be scale invariant, one needs n = — 1 or 2. The latter case with 
n > g is subject to unstable growths of curvature perturbation on large scales; hence we only 
consider the former case with n = — 1 in this paper. 

Therefore, in order to have a scale invariant power spectrum, it suffices to find functions 
a(y) and c s (y) such that q = ay/ e/c s oc l/(—y). In general, the relation q oc {—y) n can be 
written as 

dlog<? . . 

dlog(-y) 

where n = —1 in the scale invariant case. Define the following flow parameters in a similar 
way, 

dloga aHy 

P= .A = — . ( 9 ) 



dlog(-y) 

d\og c s 
dlog(-y)' 

dlog(aH) <ilog|p| 



dlog(-y) dlog(-y) 



(10) 

+ 8-1. (11) 



Note that the more commonly defined parameter s = c s /Hc s can be expressed as s = s/p. 

The parameter p measures the expansion or contraction rate of the universe. Since y 
goes from —00 to 0, p < corresponds to an expanding universe, whereas p > corresponds 
to a contracting universe. The parameter s measures the time dependence of the speed of 
sound c s . c s = 1 corresponds to s = 0; for s < the speed of sound c s increases with time, 
whereas for s > it decreases with time. The parameter r measures how the comoving 
Hubble horizon 1/aH scale with the freeze-out horizon \/\q/q"\ ~ (— y). When e or c s varies 
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with time, these two horizon scales do not coincide. Note that the parameter p also measures 
the ratio between the freeze-out horizon (— y) and the "comoving sound horizon" c s /aH. 



In terms of these flow parameters, the parameter e = |(1 + w) is given by 



-H 



e 



V 



H 2 



P 



(12) 



The null energy condition implies e > 0, which will be respected in this paper. To compute 
n, define the parameter 



V 



dlog e 
dlog(-y) 



1 

p — r 



r(r — s + 1 



(13) 



It is related to the usual slow-roll parameter rj 
it is straightforward to write n in eq. (fsl) as 



dr 

dlog(-y) 

e/He by rj = fj/p. With these parameters, 



n = p+ -(fj 



To have a constant re, equations (11) and (13) imply the following equations 



dp 



dlog(-y) 
dr 

rflog(-y) 



-p(r - s + 1), 

-p(2p — s — 2n) + 2pr — r(r + 2n + 1). 



(14) 

(15) 
(16) 



On the left hand side, the derivative is with respect to the "log time" — log(— y), which runs 
from — oo to oo. Given the value of s, these two equations determine the flow in the parameter 
space (p, r, s) . 

For example, in the case where c s = const, the parameter s = 0. In the pr-plane with 
s = 0, the flow lines according to eqs. (15, 16) for n = —1 are shown in figure [T] The shaded 
region corresponds to e = (p — r)/p < 0, which is forbidden by the null energy condition. 
There are three fixed points at (p,r) = (—1,-1), (0,0), and (0,1) respectively. The one at 
(0, 1) behaves as a sink where nearby trajectories converge. The one at (0, 0) is a saddle 
point where trajectories approach from both sides but eventually move away. The other 
saddle point at (—1, —1) is attached to a separatrix which, together with the shaded region, 
divides the plane into three sectors. The flow lines in the bulk of each sector are qualitatively 
similar. 

Every flow line in the figure corresponds to a cosmological model that yields scale invariant 
perturbations. But typically the time dependence of the flow parameters along a trajectory is 
very complicated. The more interesting trajectories are those that come near the fixed points. 



There the flow velocities given by eqs. (15 16) approach zero, hence the flow parameters stay 
near those points for a sustained period in log time — log(— y). Since — log(— y) ~ log A; at 
horizon crossing, during that period a large number of modes exit the horizon with a scale 
invariant spectrum. The corresponding cosmological models can be characterized by nearly 
constant values of the flow parameters, which allows simple analytic descriptions of their time 
evolution. Indeed, the three fixed points in figure [T] precisely correspond to the three well 
known types of scale invariant models with c s = 1. Let us identify each of them. 

The fixed point at (—1,-1) corresponds to inflation [H [21 [3]. On the point (—1,-1) 
the parameter e = 0, which corresponds to the false vacuum inflation where H = const 
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Figure 1: Flow lines for n = — 1 in the pr-plane with s = 0. The three solid circles are 
the fixed points representing inflation (—1,-1), adiabatic ekpyrosis (0,0), and apex (0,1) 
respectively; the dashed line is a separatrix. Each line segment between adjacent arrowheads 
on the flow lines represents 0.1 log time interval, or 0.1 e-fold of scale invariant modes. 

until the universe tunnels to the true vacuum. In that case aH = l/(— r), where r is the 
conformal time that equals the sound horizon time y for c s = 1, hence p = aHr = — 1 and 
r = clHt(1 — e) = —1. In slow-roll inflation, however, e is small but nonzero; the energy of 
the inflaton field <p is dominated by a nearly flat potential V((f>). The slow-roll approximation 
implies « — V^/SH rj {—V^/V)-\JV /2>, and e « \{^f^) 2 "C 1 which is assumed to vary 
only slowly. Therefore 

*' = J H *~J&<3&m~Wh (17) 

hence a ~ e^~ V ' V '^ and 

ydt r e WWd4> e(W / e (v/ W (z^^~_ J_ + er fl8 ) 
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where an integration by parts is performed and the integration constant is chosen to be zero. 
Consequently, p = aHr — 1 — e and r = aHr(l — e) « — 1 + 0(e 2 ), indeed close to the 
fixed point. Inflation ends when the potential becomes steep and the field speeds up, which 
lead the flow parameters to move up across the r = (e = 1) line and stop producing scale 
invariant modes. 

The fixed point at (0,0) corresponds to the adiabatic ekpyrotic model [HE]. This model 
relies on a slow contraction phase with H ~ Hq + c/t, where Hq < and c<Cl. During the 
period between t^ eg m 1/Hq and t en d ~ c/Hq, H rj Hq but H ~ —1/t 2 , hence e increases 
rapidly as ~ 1/t 2 . Also within this one Hubble time a ~ const, hence the conformal time 
r ~ t. Therefore q 2 ~ e ~ 1/t 2 , which guarantees a scale invariant power spectrum. The 
number of scale invariant modes is given by N m log |ibegAend| ~ log(l/c). To see that during 
this period the flow parameters stay near the saddle point, note that p = aHr « Hot, and 
p — r = aHre —tH/H « c/Hot. Hence for ibeg C f <C i C nd we find < p <C 1 and \r\ <C 1. 
Note that (p — r) corresponds to a growing eigenmode of eqs. (15, 16) near the saddle point 
(0,0); it starts extremely small at ibeg = 1/Hq, (p — r)b eg ~ c ~ e~ N <C 1, and becomes 
(p ~ r )end ~ C(l) a * *end ~ c/Hq. We have shown that the adiabatic ekpyrotic mechanism 
takes place near the point (0, 0) on the p > side. Notably, this fixed point may also be 
approached from the p < side, which corresponds to a slow expansion [9j [18] that produces 
scale invariant perturations in a similar way. 

Finally, the fixed point at (0,1) corresponds to the apex model, first discussed in [9]. 
In this model the universe slowly transitions from expansion to contraction. But unlike 
in [9], here we describe the apex model as a distinct mechanism that is independent from 
the adiabatic ekpyrotic slow expansion. When the universe gradually stops expansion to 
reach an apex, the scale factor a slowly comes to a halt, and the Hubble parameter H has 
to drop from positive to zero. Therefore, near the apex we may assume that H ~ Hot, 
where the constant Hq < and t goes from negative to 0. During the last Hubble time 



when t > ibeg = — l/v — Hq, the scale factor a ~ const, hence r ~ t. The scale invariant 
perturbations are produced by q 2 ~ e ~ —Hq/H 2 ~ 1/t 2 . The flow parameters are given 
by p = aHr w Hot 2 , and r = aHr(l — e) 1 + Hot 2 . They approach the sink point (0, 1) 
asymptotically as t — > 0~ . It seems that, as the log time — log(— r) — > oo, an infinite number 
of scale invariant modes would be produced. Nevertheless, before reaching the point (0, 1) at 
which e — > oo, our analysis based on the quadratic action ([T]) should already be expected to 
break down due to strong coupling [10J. Therefore the number of scale invariant modes should 
be finite. It is also important to have an asymmetry in the cosmic evolution before and after 
the apex, so that the scale invariant modes are not undone in the subsequent contraction 
phase. 

We have represented all scale invariant cosmological models with c s = 1 by the flow 
lines on the s = plane in the parameter space. They correspond to solutions a(r) of the 
generalized Emden-Fowler equation imposed by the requirement of a scale invariant power 
spectrum [ID] . 

-Yoci. (19) 
a J t z 

Our diagrammatic representation provides a clear way of classifying models by the flow lines 
and searching for simple analytic solutions near the fixed points. The main advantage of our 
approach will become clear when we consider models with a time varying speed of sound. In 
that case the requirement of a scale invariant power spectrum alone does not suffice to fully 
determine the time dependence of both functions a(y) and c s (y), hence it is not possible to 
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Figure 2: Positions of the four fixed points in scale invariant models (n = —1), projected on 
the pr-plane for different values of 3: the adiabatic ekpyrosis (solid blue circle) at (0,0) for 
all 3, the decelerated expansion (red circle) at (0,1) for all 3, the inflation/deflation point 
(thick orange line) at (3—1,3—1), and the tachyacoustic expansion (dashed green line) at 
(| - 1,5- 1) for < 3 < 2. 



find solutions without imposing further constraints. On the other hand, in terms of the flow 
parameters, eqs. (15, 16) still hold except that 3 is now another parameter that may vary. 
To look for new types of scale invariant cosmological models, we shall generalize our analysis 
to the whole parameter space (p, r, 3) . 



3 Fixed points 

In the general case, the cosmological models will be represented by flow lines in the parameter 
space (p, r, 3). Since the types of cosmological models are largely determined by the positions 
and properties of the fixed points, we will first locate the fixed points and then analyze the 
flow lines around them. 



Despite the unknown dynamics of 3, the fixed points of the flow equations (15 16) must 
satisfy 

f -p(r -3+ 1) = 0, 

{ (20) 
I —p(2p — 2n — 3) + 2pr — r(r + 2n + 1) = 0. 

For each value of 3, there are four solutions to these equations, namely (p, r) = (0, 0), (0, —2n— 
1), (3 — 1,3—1), and (3/2 + n,s — 1). Each solution draws a fixed point in the parameter 
space, which traces out a line as 3 varies; their positions projected on the p r-plane are shown 
in figure [2] for n = — 1. The first three fixed points are generalizations of the ones in the c s = 1 
(3 = 0) case to models with time varying c s (3^0). They will be referred to as the adiabatic 
ekpyrosis point (0,0), the decelerated expansion point (0, 1), and the inflation/ deflation point 
(3—1, 3—1). The fourth fixed point (3/2 — 1, 3 — 1) corresponds to the tachyacoustic expansion 
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scenario [HI [121 113 1H]- It is degenerate with inflation at s = 0, but otherwise is a distinct 
scenario that can generate scale invariant curvature perturbations in an expanding universe 
with a decreasing speed of sound. Incidentally, for n = 2 the fourth fixed point (s/2 + 2, s — 1) 
corresponds to a contracting universe with a growing speed of sound described in |16j . 

To learn the properties of the fixed points, we shall analyze the flow lines near those points, 
as in the s = case. In the more general case, the parameter s may be nonzero and even time 
dependent. However, in the larger parameter space (p, r, s), the dynamics of the parameter s 



is not constrained by the requirement of a scale invariant power spectrum that led to eqs. ( 15 



16). For a practical model such as a scalar field with a particular Lagrangian, the speed of 
sound c s is determined by the field configuration. Hence generally there is a relation between 
the parameters s and p, r, as shown by a few examples in section [4} Here, without being 
constrained by the underlying model, we will phenomenologically explore different situations 
in which s can vary with respect to p and r. 

In particular, consider flow lines on a surface s(p, r) in the parameter space that goes 
through one fixed point. Then near that point a small deviation in s is given by 

ds ds 

5s « — bp + — 5r = a 5p + B 5r, (21) 
op or 



where (a, B) is the gradient of the function s(p, r) at the fixed point. The flow equations ( 15 



16) can be expanded to linear order in 5p, 5r near that point, setting n = — 1 
d5p 



dlog(-y) 
d5r 



{-ap + r - s + l)6p- (1 - B)p5r, (22) 
((4 - a)p - 2r - s + 2) 5p - (-(2 + B)p + 2r - 1) 6r. (23) 



dlog(-y) 

The eigenvalues to this linear system are given by 
j Ai + A 2 = (2 + a + B)p - 3r + s, 

| AiA 2 = (3a + 4/3 - 4)p 2 + (3s + a + 8 - 4)p - (2a + 38)pr + (2r - l)(r - S + 1). 

These eigenvalues determine whether a fixed point is a source, a sink, or a saddle point. 

In the following we analyze each of the four fixed points and derive approximate solutions 
of the nearby flow lines. The corresponding types of cosmological models are described and 
verified to create scale invariant curvature perturbations. These models rely on having a 
manifestly time dependent speed of sound; constraints from the strong coupling problem and 
superluminality are discussed in section [5} 

3.1 adiabatic ekpyrosis: (p, r) = (0,0) 

This fixed point stays in the same position for any value of s. Like in the s = case, 
the corresponding cosmological model is either a slow contraction if < p <C 1 or a slow 



expansion if —1 <C p < 0. The eigenvalues from eqs. (24) can be explicitly found to be 
Ai = 1, A 2 = s — 1. Because the first eigenvalue is positive, this fixed point can be either a 
saddle point or a source of flow depending on the value of the second eigenvalue. Specifically, 
it is a saddle point if s < 1, which is qualitatively the same as in the case s = discussed in 
section [2} However, for s > 1 we find a new situation where this point behaves like a source. 
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In this new situation with 3 > 1, the flow lines emerge from the point (0,0) in the 
asymptotic past y —> — oo. In that limit, since p, r — > 0, both a and H approach constant 
values, a r& (2_oo stud H ~ H—^. Take ot — j3 — for simplicity, so that c s ~ ( — y) s • Then in 
the asymptotic past, 

t= Ji du ~Li^ dy ~F^> (25) 

which goes from to oo. That is, the asymptotic past corresponds to finite physical time t. 
To see the cosmic evolution in this situation, it is better to solve the flow equations and 



express a and H in terms oft. Near (p,r) = (0,0), the eigenmodes of eqs. (22 23) are given 



bypocl/(— y) s 1 and r — p oc l/(— y). Using the definitions of p and r, one finds 

log (^)~H^~ ( ' (26) 

Therefore, in the limit t — > 0, 



a«a_ 00 (l + F_ 00 t), (28) 

H « F-oo (l - C(fl-_ 0O t) 1/(S - 1) ) , (29) 

^-(-yf-i-^ 5 - 1 ), (30) 

where we assume C ~ O(l) > 0. For t <C l/|S r _ 00 | it is approximately true that a ~ a-oo 
and H ~ i?-oo- Under this approximation, 

" 7/ ( * (F-oot)^^*- 1 ) ~ (-y) 5 " 2 , (31) 



which guarantees the condition q 2 ~ e/c s ~ l/(—y) 2 for generating scale invariant curvature 
perturbations. 

A particularly interesting example is when s = 2. This corresponds to a cosmological 
model in which H r* -ff-oo — CH^^t. For < £ <C l/|i^_oo|, it gives a « a_oo and 
e w C = const. If the speed of sound satisfies c s ~ 1/t 2 , then one finds (— y) ~ 1/t, hence 
q 2 ~ l/c s ~ l/(— y) 2 . That is, the scale invariance of the power spectrum is achieved by 
having a uniquely time dependent speed of sound. This model can be considered as dual 
to either inflation or adiabatic ekpyrosis in which the l/(— y) 2 dependence of q 2 = a 2 e/c s 
comes solely from either a or e. Like in adiabatic ekpyrosis, this mechanism works in either 
a contracting or an expanding universe according to whether -£f_oo is negative or positive. 

3.2 decelerated expansion: (p, r) = (0, 1) 

This fixed point represents a decelerated cosmic expansion, since near this point p < and 



e = {p — r)/p S> 1. The eigenvalues from eqs. (24) are Ai = —1, A2 = 3 — 2, which again do 



not depend on a, f3. For s < 2 this fixed point is a sink, the same as the apex model in the 
s = case. But for 3 > 2 it becomes a saddle point, which gives rise to a new scenario. 

In this new scenario with 3 > 2, the flow lines do not asymptote towards p = but move 
away from it instead. Hence the universe does not approach an apex. Indeed, near the point 
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(0, 1), a w const and H ~ (—?/); but unlike in the apex model where H — > as y — > , here 
the flow line passes by the saddle point while y remains finite. The unstable eigenmode of 



eqs. (22 23) is given by p oc — 1/(— y) s ~ 2 , which implies e ~ — 1/p oc (—y) s ~ 2 - For s > 2, e 
decreases with time rather than increases. 

To see the cosmic evolution in physical time t, assume again a = (3 = so that c s ~ (—y) s . 



Then for a ~ const one obtains t ~ \/{—y) s 1 as in eq. (25), where t > and increases. Thus 



the cosmological model can be described by a decelerated expansion with 

ff-r 1 ^" 1 ', c s ~ t-W-i). (32) 

We can check that if ~ —^-s/O- 1 ) anc j hence e ~ t-( s - 2 )/( s - 1 ) (—y) 5 - 2 ^ ensuring g 2 



e/c s ~ l/(—y) . To be consistent, since by eq. (32), 



logf-W*- 2 )/^- 1 ), (33) 
\a y 

the approximation a ~ const holds for sufficiently small t. This approximation breaks down 
at finite time before H reaches 0, showing again that it is different from an apex model. 

There is a marginal case where s = 2 + 2po with — 1 <C po < 0. In that case p w po = const, 
hence a ~ (—y) po - Since c s ~ (— y) 2+2po , the physical time t ~ l/(— y) 1+po > 0. The cosmic 
evolution is given by 

This is a slow expansion with constant e = — (1 +Po)/Po 3> 1. Therefore g 2 ~ a 2 /c s ~ 
l/(— y) 2 , which produces scale invariant curvature perturbations. This model describes an 
extremely decelerated expansion with a rapidly decreasing speed of sound, which is the £>1 
limit of the tachyacoustic model to be discussed in section |3.4| Interestingly, this model has 
a contraction analog with constant e>l, represented by a flow line in the p > half plane. 



This new "acoustic ekpyrotic" model will be discussed in section 4.1 



3.3 inflation/deflation: (p,r) — (s — 1, s — 1) 

This fixed point describes a de Sitter universe with e = (p — r)/p = 0. The position of this 
fixed point depends on s. For s < 1 it lies on the border p = r < in the left half plane, 
the same as inflation; but for s > 1 it moves to the right half plane with p = r > (the 



case s = coincides with adiabatic ekpyrosis). The eigenvalues from eqs. (24) are Ai = s, 
A2 = (a+/3— — 1). Depending on the value of a + /3, this point can be a saddle, a source, 
or a sink, as listed in table [T] 

Here we consider the new situation when s > 1 and p = r > 0, which corresponds to 
an extremely rapid cosmic contraction, or "deflation". Specifically, near this fixed point, 
p ~ s — 1 and r — p~ 0, hence a ~ (— y) s_1 and if ~ const. The eigenmode that corresponds 
to the first eigenvalue is (p — r) oc l/{—y) s . Therefore e ~ (p — r)/(s — 1) ~ l/(— i/) s , which 
is crucial to ensure q 2 = a 2 e/c s ~ l/(— y) 2 . Now that e is rapidly increasing, it is clear that 
the flow parameters can only stay near the deflation point for a finite time. 

For a ~ (— y) s ~ l and c s ~ {—y) s , the physical time t is given by 
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i.e. proportional to the number of scale invariant modes. In physical time t, if H ~ Hq < 0, 
then a ~ e Hot . Comparing that to a ~ (— y) s_1 , one finds 

t«^log(-y), or (-y)^e^- 1 ). (36) 
Hq 

The crucial part of this model is to have a rapidly growing e and a rapidly decreasing c s , 

e «^° raeoe -A"ot c s ~ (-y)* ~ e sV^. (37) 

(-vr 

The de Sitter approximation e <C 1 is valid until t en( j ~ ^(-Hq) ^ we assume that 

eo -C 1, then there can be many numbers of Hubble times before t e nd- 

Therefore, this deflation model is described by a rapid cosmic contraction with 

H^H (l + ^e e-A Hot ), (38) 



which satisfies H Hq for t <C t e nd- The scale invariance of the power spectrum is achieved 
by an interplay between 

a~e Hot , e^e e-^ Hot , c s ~ e ^ Hot . (39) 

So far the dependence of the model on a, f3 has been neglected as we implicitly assumed 
a nearly constant s. For a + f3 > 1 the fixed point is actually a source where the flow line 
emerges in the asymptotic past, so the above description works even better at early times 
t —7- — oo. However, for a + f3 < 1 the fixed point is a saddle point instead. In that case, when 
extrapolating backwards in time, the flow parameters would deviate from the saddle point at 
certain time tbeg- This time is controlled by the second eigenvalue, ibeg ~ (i- a -p)H 
The model described above is valid between £beg and i en d to produce roughly N ~ log(^) 
scale invariant modes. 



Note that eqs. (38, 39) also apply for s < 1 and Hq > 0, in which case they describe an 
exponentially expanding universe with a time varying speed of sound, or "acoustic inflation" . 
This model is less interesting since inflation with e <C 1 and a constant c s can already 
produce scale invariant perturbations. In contrast, the deflation model relies on having a 
rapidly decreasing c s as well as a small but increasing e. 

3.4 tachyacoustic expansion: (p,r) = (s/2 — 1, s — 1) 

The position of this fixed point also depends on s. Moreover, since e = (p — r)/p = s/{2 — s), 
the null energy condition e > requires < s < 2 (the cases 5 = 0,2 coincide with inflation 
and apex respectively). Thus in the pr-plane, this fixed point lies on the line segment 
r = 2p + 1 for — 1 < p < 0, as shown in figure [2] Since p < and s > 0, it describes an 
expanding universe with a decreasing speed of sound. 

Here we show that this model precisely corresponds to the tachyacoustic expansion sce- 
nario [HI H21 [131 E]. Specifically, at the point (p, r) = (s/2 — 1, s — 1) for certain value of s, 



the model is described by a ~ (—y) s ^ 21 , H ~ {—y) s ^ 2 , and c s ~ (—y) s . From eq. (14) one 
finds fj = 0, consistent with e = 5/(2 — s) = const. Therefore q 2 ~ a 2 /c s ~ l/(—y) 2 , ensuring 
a scale invariant power spectrum. 
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In this model, the physical time t is given by 

dy 



t 



-dy 



(40) 



_ y )S/2+l {-yf/ 2 ' 
Thus in physical time t, a ~ 

t (2-s)/5 9 H = (2- 5)/(5t), and c s ~ 1/t 2 . We can check that 
e = —H/H 2 = 5/(2 — 5) as expected. In terms of the parameter s = c s /Hc s = s/p, this 
relation can be expressed as s = — 2e, which agrees perfectly with the tachyacoustic scenario 
\12~[ [T3"l [T4l [T6] . This scenario is emphatically different from inflation since here e need 
not be small at all. In particular, the limit e> 1 corresponds precisely to the marginal case 
in section I 



For flow lines near this fixed point, the eigenvalues in eq. (24) become 



Ai + A 2 = («+£_i)S-(a + 0-l), 
AiA 2 = i(f + /3-l)5(2-5). 



(41) 



Therefore this point is a saddle point if § + /3 < 1, otherwise a source or a sink depending 
on whether - l)5 - (a + f3 - 1) > or < 0, as listed in table [l] 



fixed adia ekpy 
points (0, 0) 



source 



s > 1 



saddle 



s < 1 



sink 



decel exp 
(0,1) 



~s>2 



s<2 



inf / def 
(5-1,5-1) 



0<5<l,a + /3<l 
5 > I, a + P > 1 



5<0,a + /3<l 



0<5<l,a + /3>l 



s>l,a + /3<l 



5<0,a + /3>l 



tachyacoustic 
(|-1,5-1) 



, 2 



§+/?> 1, 

1)5 > (a + /3-l) 



f + /3 < 1 



f +/3> 1, 



^-1)5< (a + /3-l) 



Table 1: Positions and properties of the fixed points for flow lines on a constrained surface 
s(p,r) with local gradient (a, (3). 



4 Examples 

We have analyzed the flow lines near the fixed points under certain constraint function s(p, r) 
with local gradient (a,/3). For concreteness, consider examples that are motivated by scalar 
field models. Suppose the universe is dominated by a scalar field (j) with Lagrangian C = 
P(c/),X), where X = —^(dcj)) 2 . It is shown in appendix [X] that, for a given function P(4>, X), 
the parameter r/ = e/He can be locally determined as a function of e and s. Such a function 
7/(e, s) suffices to provide the constraint between the parameters p, r and 5. 

4.1 example I: rj = ms 

Consider first the A;-essence like model in which the Lagrangian is factorizable as C = 
F((p)P(X) |19j . In this case both the equation of state w and the speed of sound c s de- 



13 



pend only on X , namely 



P 



w 



2XP 



P 



2XP XX + P x 



(42) 



Therefore the parameter e = 1(1 + w) is locally a function of c s , which is determined by the 
form of the function P(X). 

As a simple example, take the function e(c s ) to be a power law, e oc c s m , with a constant 
exponent m. Such a relation has been considered in different models. The case with e oc (? a 
(m = 2) appears in nonsingular bouncing models where a ghost condensate field violates the 
null energy condition. The case with e oc c s (m = 1) is considered in |15j for which some 
contribution to the non-Gaussianity is scale invariant, so as to avoid the strong coupling 
problem. And the case with e = const (m = 0) leads to the tachyacoustic model. 

This above relation e oc c s m is equivalent to having rj = ms, a special case in which r\ 
only depends on s. In terms of the flow parameters, it can be expressed as fj = ms. Then 



eq. (14) implies, for m/1, 



2n 



s=- p + -. (43) 

1 — m m — 1 

For a constant value of n, this function s(p, r) determines a plane in the (p, r, s) parameter 
space that is parallel to the r-axis and goes through the point (p, s) = (n, 0). Comparing it 
to eq. (plL we identify a = /3 = 0. 

To find the fixed points and their properties in this example, we follow the analysis in 
section [3] and take n = —1. The adiabatic ekpyrotic fixed point is at (p,r,s) = (0,0, jz^), 
which is a source if s = > 1, i.e. — 1 < m < 1, or a saddle point otherwise. The 

> 2, i.e. 



1— m 



decelerated expansion point is at (0, 1, tz^)i which is a saddle point if s 
< m < 1, or a sink otherwise. The inflation point is at (—1,-1,0) regardless of m; since 



1— m 



0, it is a sink if a + ft 



> 1, i.e. 



i_ m . 1 < m < 1, or a saddle point otherwise. 

For m 7^ 0, the tachyacoustic point is also at (—1,-1,0), degenerate with inflation; but for 
0, it can be any point (s/2 — 1, s — 1, s) with < s < 2 — this is just the tachyacoustic 



in 



model discussed in section |3.4| These fixed points are summarized in table [2] 



m 

(a = 5^/3 = 0) 


adia ekpy 
(0.0. l-m) 


decel exp 

(0.1. i-J 


inflation 
(-1,-1,0) 


tachyacoustic 
(| - 1,5 - l,s),0 < s < 2 


m < —1 


saddle 




saddle 


/ 


-1 < m < 




sink 






m = 


source 




sink 


sink 


< m < 1 




saddle 




m = 1 


/ 


/ 




/ 


m > 1 


saddle 


sink 


saddle 





Table 2: The fixed points in the model e oc c s m (rj = ms). 



To be complete, the m 
finds p - 



1 case is treated separately here. Since fj = s, from eq. (14) one 
1, then eq. ( 15 ) implies r = s — 1. The null energy condition e = 1 + r > further 
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2 



1 



r 



-1 



-2 




-2-10 1 2 

P 



Figure 3: Flow lines for models with e = const. The adiabatic ekpyrosis point forms a source 
at the origin, and the tachyacoustic expansion points form an attractor line connecting the 
decelerated expansion point and the inflation point. In this and the following figures each 
line segment between adjacent arrowheads represents 0.5 e-fold of scale invariant modes. 



restricts the parameter space to a half line {p = — 1, r = s — 1 > — 1}. In this restricted 
parameter space, neither the adiabatic ekpyrosis nor the decelerated expansion exists, whereas 
inflation can still happen at (p, r, s) = (—1, —1, 0) (degenerate with the tachyacoustic point). 
This point is a sink which the flow parameters quickly approach, as can be found by solving 



eq. (16). This model is described by a oc l/(— y) and c s oc e = —l/log(y/yo), where yo < y < 
0, as found in [T5] . 

Here we present in detail the case with m = 0. The parameter space is constrained to the 
plane s = 2p + 2. The flow lines have radial trajectories, as shown in figure [3] It is clear that 
the adiabatic ekpyrosis point at the origin is a source, whereas the tachyacoustic points form 
an attractor line segment, with the decelerated expansion and the inflation points at the two 
ends. Since the flow lines in the left half plane quickly approach the attractor points, the 
corresponding cosmological models can be effectively described by the tachyacoustic expan- 
sion in section [3T4| On the other hand, the flow lines emerging from the origin in the right 
half plane correspond to contraction models with a decreasing speed of sound. In particular, 
the trajectories near the lower half r-axis have e = const S> 1 — they represent precisely the 



acoustic ekpyrotic contraction mentioned in section 3.2 

To describe this "acoustic ekpyrotic" model, we analytically solve the condition q 2 
a 2 e/c s oc l/(— y) 2 with the constraint 



a 2 



const. (44) 
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The result is 



a{y) oc 



1 



i 

e+1 



(45) 



Hfc) H) 

where y < y c < 0; it describes a contracting universe that crunches at y c . Accordingly, the 
flow parameters can be found as 

d log all 



P 



dlog(-y) e + l(JL)-r 



r = (1 — e)p 



1 



1 



e + 1 (^)-l' 



(46) 
(47) 



(Incidentally, y < < y c would describe the flow lines between the origin and the attractor 
line, and y c < y < the flow lines to the left of the attractor line.) Here these parameters 
vary with time as they flow along the trajectories, unlike near the fixed points where they 
lead to simple power- laws. 

Since c s oc a 2 (—y) 2 , the physical time t is given by 



t oc 



ady 
o 2 (-y) : 



oc 



1 



1 



Hfc) H) 



(48) 



which goes from to oc —(—y c ) e ^ e+1 ' to as y goes from -co to y c . From this we recover 

aoc(-t) 1 / £ , H=—, (49) 



which exactly describes an ekpyrotic contraction for e> 1. The scale invariance of the power 
spectrum relies on a decreasing speed of sound, 



c, oc 



(-i)2A 



1 



.to) 



1 



2 " 



(50) 



where the approximation is due to a ~ const during the ekpyrotic phase. Note that as y — > 
-co, or t — > to, both a and -ff are finite. Near that limit the flow parameters asymptotically 



emerge from the adiabatic ekpyrotic point, which is the situation described in section 3.1 
hence in that limit the exact solutions above can be approximated by the results there. 

4.2 example II: rj = ae + b 

Another example is motivated by a particular eft dependence of the function F((f>) in the k- 
essence model, namely F((f>) oc 1/0 2 [20l[l9]. In this case, it is shown in appendix [a] that the 
parameter r\ can be determined as a function of e only. 

For simplicity, take 77(e) to be a linear function, 77 = ce + b, with some constants b, c. 
Since by definition rj = e/He = dloge/dloga, such a function 77(e) can be used to study 
the running of the parameter e. Namely, in an expanding universe if c < then e should 
approach the fixed point e* = —b/c, whereas in a contracting universe the fixed point is an 
attractor if c > 0. 



Using eqs. (12) and (14), the relation 77 = ce + b can be expressed as 



+ c + 2) p 



cr 



2n, 



(51) 
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Figure 4: Flow lines for c < < b (here c = —l,b = 1). The adiabatic ekpyrosis point 
forms a source at (0,0), whereas the decelerated expansion point (0,1) and the inflation 
point (—3,-3) are two saddle points connected by a separatrix. The tachyacoustic point 
(— i, 0) on the separatrix forms a sink where the flow lines converge and e approaches e* = 1. 



which also determines a plane in the parameter space (p, r, s) for constant values of n. In 
comparison with (21 ) we identify a = b + c + 2, f3 = —c. 

Following the general analysis, we look for fixed points with n = — 1. The adiabatic 
ekpyrotic point is at (0, 0, 2), which is a source since s > 1. The decelerated expansion point 
is at (0, 1, 2 — c), which is a sink if s = 2 — c < 2, i.e. c > 0, or a saddle point otherwise. The 
inflation/deflation point is at (jEct, j+x ; 5+j) for b 7^ — 1, which is a sink if s = 5^ < an d 
a + /3 = 6 + 2>l, i.e. — 1 < b < 0, or a saddle point otherwise. The tachyacoustic point 
2b ) provided < s = < 2; it is a saddle point if b < < c, or a sink if 



is at (b^ c i gzf , fc _, 



c < < 6. These cases are listed in table [3j 



a = b+c+2 
P = -c 


adia ekpy 
(0,0,2) 


decel exp 
(0,1,2-c) 


inf / def 
V6+l» 6+1 ' 6+1 J 


tachyacoustic 

/ c 6+c 2b \ 
Vb-c' b-c' b-c/ 


source 


V6,c 


/ 


/ 


/ 


saddle 


/ 


c < 


6<-lor6>0 


b < < c 


c = 


6 < 


6 = 


c> 


sink 


/ 


6 > 


c < 


c < < 6 


c> 


-1< 6 < 



Table 3: The fixed points in the model rj = ce + b. 
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Figure 5: Flow lines for — 1<6<0<c (here c = 1,6 = — g). Compared to figure [ij the 
adiabatic ekpyrosis point (0, 0) is still a source; both the decelerated expansion point (0, l) 
and the inflation point (— |, — |) become sinks, whereas the tachyacoustic point (— |, — ^) is 
no longer an attractor. The flow lines in the right half plane will approach the slope 1 — e* = | 
further away from the origin. 



As an example, consider the case c < < b. The flow lines projected on the pr-plane is 
shown in figure |4 There is a source point at (0, 0), two saddle points at (0, 1) and (gq^j, 
connected by a separatrix that goes through a sink point at (jz^, 533 )■ In the right half plane, 
all trajectories emerge from the origin and flow away to infinity; the cosmic contraction near 



that point can be well described by the generalized adiabatic ekpyrotic model in section 3.1 



The left half plane is divided by the separatrix into two sectors. Trajectories in the right 
sector emerge from the origin and flow towards the sink, whereas trajectories to the left of the 
separatrix come from infinity and also flow to the same point. Hence in an expanding universe 
e always approaches its value at the sink point, e = (p — r )/p\ sink = —b/c = e*, as predicted 
by the running according to rj(e). The cosmic expansion near that point is described by the 



tachyacoustic model discussed in section 3.4 



Similarly, an example for the case — l<6<0<cis shown in figure [5} Here in the left 
half plane, the tachyacoustic point becomes a saddle point while the two endpoints of the 
separatrix become sinks instead, hence e is repelled away from the value e* in an expanding 
universe. Meanwhile, in the right half plane, the trajectories emerge from the origin and 
eventually go parallel to the line with a slope r/p = (b + c)/c = 1 — e*, proving that e 
approaches e* in a contracting universe. 

The case with b < —1,0 < cis even more different, as show in figure [6j In this case 
the inflation point disappears from the left half plane; it is replaced by a saddle point on 
the border of e = in the right half plane, attached with a separatrix that divides the 
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Figure 6: Flow lines for b < —1,0 < c (here c = 1,6 = —2). Compared to figure [5j the 
inflation point in the left half plane is replaced by the deflation point in the right half plane, 
which is a saddle point at (1,1). The flow lines in the right half plane approach the slope 
1 — e* = — 1 away from the origin. 



right half plane into two sectors. This is the new situation discussed in section 3.3 — near 
that saddle point, the cosmic contraction can be described by the deflation model with a 
decreasing speed of sound. Further away from that saddle point, the flow lines tend to 
approach the slope r/p — > 1 — e* just like in the case above; the cosmic contraction with an 
asymptotically constant e = e* can be described in the same way as the acoustic ekpyrotic 
model in section HUl 



5 Physical constraints 

In the simplistic analysis above, each model can potentially produce a large number of exactly 
scale invariant modes. In practice, some physical constraints should be considered that could 
affect the stability of the models and limit the range of scale invariant modes |2H [TO] . 

First of all, the period during which the scale invariant modes are created should be stable 
against the growth of spatial curvature and anisotropy. The effective densities of these two 
components scale as 1/a 2 and 1/a 6 respectively [22], whereas the dominant energy density of 
the universe scales as l/a 2e . Therefore an expansion phase with e < 1 can efficiently dilute 
away spatial curvature and anisotropy, leading to a flat, homogeneous and isotropic universe 
as observed today. Hence accelerated expansion models like inflation (with e <C 1) and 
tachyacoustic expansion (with e < 1) can solve the horizon problem as well as provide a scale 
invariant power spectrum; but decelerated expansion models including adiabatic ekpyrotic 
expansion and apex (with e > 1) must be supplemented by additional mechanisms, e.g. 
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an ekpyrotic contraction phase, to suppress inhomogeneity and anisotropy |18| . Similarly, a 
slow contraction phase with e > 3, such as adiabatic ekpyrosis and acoustic ekpyrosis, can 
automatically beat curvature and anisotropy; but a rapid contraction such as deflation (e <C 1) 
would require an extra stabilization phase, just like in the matter dominated contraction 
(e = |) model [23]. 

Besides, for models where the comoving Hubble horizon 1/aH does not coincide with the 
freeze-out horizon, the scale invariant modes that have left the freeze-out horizon must also 
exit the Hubble horizon before the standard expansion phase. In models including adiabatic 
ekpyrosis, decelerated expansion, and deflation, this can be done by having a subsequent 
ekpyrotic phase to push the modes outside the Hubble horizon, such as in [H [TB] . 

Moreover, since the cosmological models assume classical general relativity, for consistency 
the energy density of the universe should remain sub-Planckian, H 2 /Mp^ < 1. To estimate 
the size of the Hubble parameter when the scale invariant modes are generated, taking n = — 1 
in eq. Q to fix the coefficients in eq. ^ , one finds that the power spectrum ([6]) becomes 

1 _ 1 

8vr 2 M|^V - &^m^ e a 2 y 2 • 



From the measured amplitude of the power spectrum, ~ 2.4 x 10 , the sub-Planckian 
energy density criterion can be expressed as 



H 2 8vr 2 ea 2 H 2 y 2 p 2 ec 



s 



Mpj 4 x 10 8 c s 5 x 10 6 



< 1- (53) 



For models with time varying e and c s , it is also important to avoid strong coupling in 
generating scale invariant perturbations. A simple estimate and often reliable criterion is 
that the magnitude of the third order action for S3 ~ f 0(C 3 ), should be smaller than the 
quadratic action S2 in |24l [TO] . For a rapidly growing e or a rapidly decreasing c s , the 
dominant contribution to the ratio S3 / S2 can be estimated by the term [15] 



S3 ,e 



2 



c C^- (54) 

This term results in a non-Gaussianity that peaks at small scales where the observational 
bound on /tvl ~ S3/C $2 is weak. Demanding that S3/S2 < 1 when the scale invariant modes 
crosses the horizon, the weak coupling criterion can be written as 



e 2 „ 1 1 



< 



c 2 ~ C 5 x 10- 5 ' 



(55) 



The above two constraints can quantitatively limit the number of scale invariant modes 
generated in certain types of models. For example, in the tachyacoustic model with a constant 



e ~ O(l), using p = s/2 — 1 = — l/(e + 1) from section 3.4, eq. (53) becomes c s e/(e + l) 2 < 



5 x 10 6 . Combining that with eq. (54), one obtains 



7xl0- 3 6<c 5 <5xl0 6 (e + 1)2 . (56) 



e 



Since c s ~ (—y) s ~ k 2e /( e+1 ) at horizon crossing, the number of scale invariant modes is 
bounded by [15] 
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Similarly, in the decelerated expansion model described in section 3.2 using e ~ — 
eqs. (53) and (54) imply 7 x 10~ 3 < c s /e < 5 x 10 6 . Since in this model c s /e ~ y 2 ~ 



1 /k 2 at horizon crossing, the number of scale invariant modes is bound by N < \ l°s(5 x 
10 6 /7 x 10~ 3 ) « 10, just enough to encompass the observed 10 e-folds of primordial density 
fluctuations. This result is the same as for the apex model [9]. 

Finally, models with a time varying c s often run into a superluminal speed of sound, c s > 1. 
Such a superluminal speed of sound does not necessarily violate the causal structure of space- 
time |25l 126]. Nevertheless, it is important to understand how constrained the cosmological 



models become if c s < 1. It is can be seen from eq. (56) that the tachyacoustic model re- 
quires superluminality (hence the name) for generating a wide range of scale invariant modes; 
otherwise this constraint becomes 7 x 10~ 3 e < c s < 1, allowing only N < \ log(10 3 /7e) « 3 
e-folds of scale invariant modes for e ~ O(l). The same bound also applies to the particular 
adiabatic ekpyrotic model described in section [3~T] and the acoustic ekpyrotic model described 



in section 4.1 Similarly, for the decelerated expansion model discussed above, since e> 1, 
one finds 7 x 10~ 3 < c s /e < 1, posing an even tighter bound on N. Hence these models 
all rely on having a superluminal speed of sound for generating a sufficient number of scale 
invariant modes, in the same spirit as the tachyacoustic model. 

6 Summary 

We have presented a set of flow parameters that are used to depict single field cosmological 
models. Scale invariant perturbations are produced in the model if these flow parameters 
satisfy specific dynamical equations. We have analyzed the flow of these parameters by 
identifying the fixed points and their properties. Besides existing models including inflation, 
adiabatic ekpyrosis, apex, and tachyacoustic expansion, four new scenarios have emerged 
from our analysis: 

• A generalization of the adiabatic ekpyrotic model with a time varying speed of sound, 



described in section 3.1 In a particular model, the scale invariance of the power spec- 



trum completely relies on the time dependence of c s . 

An extremely decelerated expansion with a decreasing speed of sound, described in 
section I 



An exponentially rapid contraction, or "deflation", described in section 3.3 Scale 



invariant modes are generated when e increases and c s decreases exponentially. 



An "acoustic ekpyrotic" contraction, described in section 4.1 Based on the usual 



ekpyrotic contraction, it requires a particularly time dependent speed of sound. 

These new models serve as distinctive examples in which both e and c s vary with time 
cooperatively to create scale invariant perturbations. The range of scale invariant modes can 
be estimated by considering physical constraints from sub-Planckian energy density and weak 
coupling of higher order perturbations. Unless subluminality is imposed, all these models can 
adequately account for the observed range of the scale invariant power spectrum. 

The flow parameters that we used to find the new cosmological models are convenient for 
obtaining the scale invariant power spectrum as well as analyzing the physical constraints. 
We note that in [27] a different set of flow parameters are defined by expanding the Hubble 
parameter H as a function of the field (j>. Those flow parameters can be used to effectively 
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reconstruct scalar field models, but these cosmological models do not acquire scale invariant 
perturbations. It may be interesting to relate the two sets of flow parameters in order to 
construct single field models that implement the new mechanisms for generating a scale 
invariant power spectrum. 
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A rj as a function of (e, s) 

Consider a scalar field (f> with the Lagrangian C = P(<j>,X), where X = -\{d<j)) 2 . The 
homogeneous evolution on a flat Friedmann-Robertson- Walker background is determined by 
the equations of motion 

4> = V2X, (58) 
X = -6c 2 s HX - p i(f> V2X/p,x- (59) 

Here the energy density p((f>, X) = 2XP X — P, and 

H = ±y r p/3, c 2 = P x /p iX - (60) 

where the ± sign corresponds to the universe being expanding/contracting. 
The flow parameters e, s and r/ are defined by 

e =JP> S= i^? "=H-e (61) 
Using the equations of motion, they can be explicitly written as functions of ((f), X), 

= 3XP X 
e 2XP X - P ' 

6X(P x P,xx - XP\ X + XP x Pxxx) I 6X 



{2XP XX + F x ) 2 V 2XP >X - P 

■ 2 



(62) 
(63) 



(P^Pxx-XPtx + XP^xPxxxWXPxt-Pj) X(P, XX P X ^-P X P XX ^) 
P x (2XP xx + P tX ) 2 ' Px( 2XP ,xx + Px) 



QXP tX 6(P )X + XP XX ) ± (P X P,<f, + XP tXX P^ - XP 7 xP X4> ) I 6X 



2XP X -P 2XP jXX + P jX XP jX (2XP >xx + P x ) \j2XP x -P 

(64) 

Here we argue that each model P(cp,X) determines a particular function r](e,s). Indeed, 



for a given function P((f>,X), eqs. (62 63) can be locally inverted to express (4>,X) in terms 
of (e, s), which can then be used in eq. (64) to express ijasa function of (e, s), at least locally. 
Below we consider a few examples. 
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A.l P(<f>,X) = K(X)-V(<f>) 



For scalar models with P((f>,X) = K(X) — V (</>), eqs. (62 63, 64) become 
3XK X 



2XK 



x 



K + V ' 



6X(K X K XX - XK 2 XX + XK X K XXX ) 



(2XK XX + K x ) 



1± 



6X 



6XK X V 2XK x -K + V 



(65) 
(66) 



6XK 



f] 



x 



6(K x + XKxx) 

2XK >X -K + V 2XK :XX + K >x 



1 ± 



6X 



6XK tX U 2XK >X -K + V 



(67) 



For a canonical scalar field we have K(X) = X. Then those parameters are given by 
3X 



x + v 



0, T) 



6X 
X + V 



6 + 



v rl 



6X 



x v x + v 



(68) 



This case is generically degenerate since s is trivial and (</>, X) cannot be solved from e alone. 
However, there is a special case where V is exponential, V tx e c f In that case we have 
Vfi/X = cV/X = c(3 — e)/e, hence 77 can be expressed as a function of e, 



77(e) = 2(e- 3) 1± 



'2e 



(69) 



For a constant potential V = Vo, for example, we find a linear relation 77(e) = 2e — 6. 
A.2 P(4>,X) = F(cf>)P(X) 

For another type of models with P(<f>, X) = F(<j))P(X), we have 
3XP X 



2XP X - P ' 
6X(P X P >XX 

~ XP 2 XX + XPxPxxx) 
(2XP XX + Px) 2 

6(PP X - XP 2 X + XPPxx) 
(2XP x -P)(2XPxx + Px) 



1 ± 



F 3/2 



I2XP 



P 



6XP 2 



1 ± 



3t 

F 3/2 



I2XP 



x 



P 



6XP 2 



(70) 
(71) 
(72) 



In this case e is a function of X only, and Tjzce versa, at least locally. Since the factor in the 
brackets are the same for s and 17, and the other factors are functions of X only, we can write 



V = /(e) s, 



(73) 



where the function /(e) is determined by the function P(X). For example, with p(X) = 
2XP tX - P oc X~i (7 ^ |), we find /(e) = 2, and hence 77 = 2s. 

A special case is where F(4>) oc (f)~ 2 , then we have F^/F^/ 2 = const, hence both s and 77 
depend only on X. Therefore 77 can be expressed as a function of either e or s alone. 
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